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The quasistatic limit of the velocity-gauge strong-field approximation describing the ionization 
rate of atomic or molecular systems exposed to linear polarized laser fields is derived. It is shown that 
in the low-frequency limit the ionization rate is proportional to the laser frequency, if a Coulombic 
long-range interaction is present. An expression for the corresponding proportionality coefficient is 
given. Since neither the saddle-point approximation nor the one of a small kinetic momentum is 
used in the derivation, the obtained expression represents the exact asymptotic limit. This result 
is used to propose a Coulomb correction factor. Finally, the applicability of the found asymptotic 
expression for non-vanishing laser frequencies is investigated. 

PACS numbers: 32.80.Rm, 33.80.Rv 
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I. INTRODUCTION 

Keldysh-Faisal-Reiss (KFR) theories are very popular 
to describe nonresonant multiphoton ionization of atoms 
and molecules in intense laser fields (see, e.g., [H, H, H, 0] 
and references therein). While the principle concept of 
ignoring the effect of the laser field in the initial state 
and the interaction of the ionized electron with the re- 
maining atomic system in the final state is common to all 
KFR approximations, the approaches differ in the details 
of their formulation. Whereas the length gauge was used 
in the original work of Keldysh the velocity gauge 
was invoked by Reiss 6] and Faisal [7]. Historically, 
the velocity-gauge variant of KFR theory is also known 
as strong- field approximation (SFA) Q. Although this 
terminology is not consequently adopted nowadays, this 
meaning of SFA is used in the present work that discusses 
exclusively the velocity-gauge variant of KFR. 

A natural test of KFR theories is a comparison of its 
prediction in the tunneling limit with the corresponding 
one of quasistatic theories [£J, [10|. The tunneling limit 
of the length-gauge version of KFR was considered al- 
ready by Keldysh for the explicit example of a hydrogen 
atom. In the derivation he employed, however, two ad- 
ditional approximations: the one of a small kinetic mo- 
mentum and the saddle-point method (SPM). In his de- 
tailed work [ft] about SFA theory Reiss has also used 
the SPM to obtain an approximation for the general- 
ized Bessel functions ('asymptotic approximation') that 
are required for the calculation of the transition ampli- 
tude. As in Keldysh' work he also adopted an expan- 
sion assuming a small kinetic momentum in order to de- 
rive the differential ionization rate in the tunneling limit. 
The algebraically cumbersome form of the in [|| derived 
'asymptotic approximation' has motivated the develop- 
ment of a simpler form appropriate for the tunneling 
regime. This approximation is again based on the 
SPM and is referred to as 'tunneling 1 approximation'. 
It was, however, recently criticized by Bauer 12[ who 
has shown numerically that the 'tunneling 1 approxima- 
tion' is much worse than the 'asymptotic approximation' 



in the case of large kinetic momenta. Since neither Q 
nor [TlT | contain explicit asymptotic expressions in the 
limit of vanishing laser frequency (u> — > 0) , a systematic 
study of the quasistatic limit on their basis is almost im- 
possible. Recently, the different versions of KFR were 
numerically compared with quasi-static theories for the 
hydrogen atom in [l3|. It was found that SFA signifi- 
cantly underestimates the ionization rate, especially in 
the limit to — > or for very strong fields. Since both lim- 
its can be well described with quasistatic theories, a com- 
parison of them with the corresponding limit of SFA can 
provide an insight into the reasons for such a disagree- 
ment. The derivation and analysis of such an asymptotic 
expression for the SFA ionization rate is the main mo- 
tivation of this work. As is shown below, the SFA rate 
does not converge to the tunneling result, if long-range 
Coulomb interactions are present. 

The asymptotic limit of SFA for to — > is also very in- 
teresting, because it may be used for deriving a Coulomb 
correction factor by comparing this limiting expression 
with the one of quasistatic theories. Rescaling the SFA 
rate (for oj ^ 0) in such a way that it agrees with the qua- 
sistatic limit for u — ► is supposed to correct SFA for 
the otherwise neglected long-range Coulomb interaction 
between the ionized electron and the remaining ion. Such 
a Coulomb correction factor was proposed by Becker and 
Faisal (see |lj and references therein) and is extensively 
used in their atomic and molecular SFA calculations 
Note, this correction factor is in fact very large and can 
amount to almost three orders of magnitude for atomic 
hydrogen and standard parameters of intense femtosec- 
ond lasers. Although it is emphasized in [l| that the 
low-frequency limit of SFA converges to the tunneling 
result, this is only shown for the case of short-range in- 
teractions. As is demonstrated in the present work, the 
correct asymptotic limit of SFA in the presence of long- 
range Coulomb interactions differs from the short-range 
case even qualitatively, since it is proportional to u>, but 
lu independent for short-range potentials. Therefore, the 
present work also allows to directly derive an asymptot- 
ically correct Coulomb correction factor for SFA. 
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The present paper is organized the following way. Af- 
ter a brief description of the ionization rate within SFA 
in which the basic formulas and notations are introduced 
f Sec. Ill A]) , an expression is derived in Sec. lIIBl that is nu- 
merically very convenient for the calculation of general- 
ized Bessel functions and thus the SFA in the quasistatic 
limit. In Sec. Ill Cl an exact asymptotic formula is derived 
for the generalized Bessel functions in the limit u — > 0. 
In this derivation neither the SPM nor any other approx- 
imation beyond the ones inherent to SFA are used and 
it is demonstrated that the SPM yields wrong results for 
weakly bound systems or very intense fields. In Sec. Ill Dl 
two simplifications are introduced that in contrast to the 
SPM or small-momentum approximation are universally 
justified in the limit lo — > 0. This allows to derive an 
exact analytical expression of the quasistatic limit of the 
SFA in the presence of long-range interactions that we 
name QSFA. In Sec. |TTT] the QSFA is discussed for the 
example of atomic hydrogen. After a derivation of the 
parameters specific to the considered atomic system in 
Sec. IIII Al the rate obtained in the weak-field limit is dis- 
cussed and compared to tunneling models in Sec. IIII Bl 
Based on this comparison, a Coulomb correction factor 
is derived for SFA and compared to an earlier proposed 
one. The range of validity of QSFA for standard laser 
frequencies is explored in Sec. IIII Cl where also a correc- 
tion is proposed that is explicitly given for the IS state 
of hydrogenic atoms. The findings of this work are sum- 
marized in Sec. IIVI 



II. THEORY 



defined as 



L(P) = i / e iS ^dt 



(2) 



where S(t) is given with the aid of the mechanical mo- 
mentum of the electron, ir(t) = p + {F/oj) sincut, as 



S p (t) = J dt' 



(3) 



For p = pnP the function L(p) can simply be expressed 
using the generalized Bessel functions (we use for them 
Reiss' definition which differs slightly from the one of 
Faisal) as 



L(p N p) = (-i) N e^J N ^,-z/2) 



(4) 



where £ = pwF(p • F)/lo 2 and z — U p /u>. In the 
high-frequency and low-intensity (so called multiphoton) 
regime the generalized Bessel functions can be very ef- 
ficiently calculated using an expansion over products of 
ordinary Bessel functions, 



J N (a,b)= ^ JN-2m(a)J m (b) , (5) 

rn— — oo 

where only a few terms are required to yield high accu- 
racy. 

Consider now the quasistatic limit defined by to — » 0. 
Introducing the Keldysh parameter 7 = klo/F, the (in- 
verse) field parameter r = k 3 /F, and the variables 



Ionization rate 



qN = pn/k 



and 



C = P'F, 



(6) 



In the single-active-electron approximation, we con- 
sider the direct transition of an electron from the ini- 
tial bound state ^0 to a continuum state ^f p due to the 
linear polarized laser field F(t) = Fcoswi with period 
T = 2ir/ui. The total ionization rate is given in the SFA 
by 



one finds 



W SFA = (2tt)- 2 J2 PN[^ + E h 

N>Nq 



dp\L(p N p)\ 2 \V Q (p N ,p) 



(1) 



with k = \fTEh where Eb is the binding energy of the 
initial bound state ^0 with its Fourier transform ^o- 
The number of absorbed photons N satisfies N > Nq — 
(Eb + U p )/uj where U p = F 2 / (4ui 2 ) is the electron's quiver 
(pondcromotive) energy due to the laser field. Finally, 
Pn = W2(~Nuj — £b — Up) is the momentum in the final 
state for an TV photon transition. The function L(p) is 



f = ^. "0 = ^(1 + 27"). (7) 

The condition lo — > leads to 7 — > 0, whereas parameter 
r is lo independent and thus unaffected. Since the numer- 
ical values of qpf and £ are usually of the order of one, 
both arguments and the index of the generalized Bessel 
function in (U) approach to infinity. In this case it is very 
problematic to use ([S]) for numerical calculations, since 
very many terms are required and their amplitudes are 
much larger than the final result. This can lead to large 
cancellation errors. In the next subsection we solve this 
problem by a transformation of the integral © to a form 
that is more convenient for numerical calculations. 



B. Efficient calculation in the quasistatic limit 

A very efficient way for the numerical computation of 
L(pnp) in the tunneling regime is possible by means of 
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performing the integration through the saddle points. In- 
troduction of the new complex variable u = sin cot allows 
to rewrite j2]) as 



L = <j> J-(u)du 



where 
S(u) 



e iS(u) 



(8) 



(9) 



2 7 3 



^ 2 + 2 7gj vC^ + 7 2 (i + g^) 

f(v) 



dv, (10) 



f(u) = Sign[Im(M)]Vl-u 2 . 



The closed contour C c in Eq. ([8]) encloses the branch cut 
[— 1, 1] of the functions /(u), S(u), and T(u). The path 
of integration C u in Eq. (jTTJ) specifies the path around 
the branch cut starting at v = i0 + and terminating at 
v = u. Since S(u) is a multivalued function, we have se- 
lected the branch cut along the negative imaginary axis. 
Nevertheless, function J-(u) (as well as f(u) ) is analyt- 
ical in the whole complex plane except its branch cut 
[-1,1]- 

There exist two saddle points u± of S(u) in the com- 
plex plane u defined by S'(u±) — and given explicitly 
by 



u ± = lp(-X±i) 



with 



p = y/l + q%(l-C 2 ) > 1, X = QnC/p < q N 



(12) 



(13) 



We introduce the straight contours C± that go through 
the saddle points u± and are given parametrically as 

u(x, Q±) — u± + xQ±, -co < x < co (14) 

starting at x — * — oo. The values of Q± are chosen in 
such a way that the contours C± are passing through the 
steepest descent, i.e. as 



Q± = 



2i 



S"(u±) 



(15) 



where the argument of Q± satisfies — 7r/4 < argQ± < 
tt/4. 

As \u\ — + oo, the function T{u) decays exponentially to 
for | argu| < 7r/4. This allows to transform the contour 
integral §5§ as 

L = j F{u)du - <j> F(u)du = L+-L- (16) 
c+ C- 

where the integrals L± can be calculated using (fT¥j) as 



L± = 



OO 

J Q±F{u± + Q±x)dx 



(17) 



The transformations above could also be used in the con- 
text of SPM, where the integration in (fl~7|) is performed 
in an approximate way using an expansion of S(u) at 
u = u± (see next subsection for more details). However, 
the expression obtained for L within the SPM is only 
approximate. Since our intention is to perform an ex- 
act calculation of L (within a controllable precision), the 
integration in (|17[) is done numerically using Gaussian 
quadrature. Moreover, it is sufficient to calculate only 
L + . Indeed, using 



/(«*) 

JS(u') 



= -/», Ql = Q T , 
= (-l^exp^je^")}* 



(11) one obtains 



L- = -(-l) iv exp[2^]L; 
Substituting $0$ into (JTBJ) yields 

L = L+ + (-l) N exp [2&] L* 



(18) 
(19) 

(20) 
(21) 



Introducing the absolute value C and the argument f2 of 
L + we obtain for the generalized Bessel function 



and 



J N (£,-z/2) = 2£cos{£-n- Nn/2) (22) 



\L\ 2 = 2£ 2 [1 + cos(2^ - 2f7 - TVtt)] . (23) 



We stress that no approximations have been done. The 
highly oscillatory integral has only been modified to 
a form that is much more convenient for the numerical 
calculation of \L\ 2 and will be used in the present work to 
obtain numerical values of the SFA ionization rate W^sfa 
(UJ) that serve as a reference for the QSFA derived below. 



C. Generalized Bessel functions in the quasistatic 
limit 



In order to derive (in Sec. HID]) an analytic expression 
for the SFA rate in the quasistatic limit, it is required 
to first find the exact limit of \L\ 2 for 7 — > 0. It follows 
from (fT2|) that u + — > in the limit 7 — > 0. Function f(u) 
is then nearly 1 in the interval (0, u+) and S(u + ) can be 
calculated using the Taylor expansion of / _1 (u) at u = 
for Im(u) > 0, 

/- 1 («) = i + y + ^-+- (24) 

Substitution of ([24]) into (| 1 1 1) and integration yields 

lS{u+) = — « + 0( 7 ). (25) 

Performing the Taylor expansion of S(u) at u = u + gives 

V2 



iS(u + + Q + x) — iS(u + ) ~ x +i- 



(26) 
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where 



'2/(u+) 



TP 




0(7 3 ) 



(27) 



has been used. Due to the smallness of Q+, function 
e *s{u) decays fastly in the vicinity of u + . Since the ex- 
pansions (f2"6")l and (j2"4")l are expected to be valid in this 
region, the integrand of L + can be rewritten as 



Q±T{u± + Q±x) — 7 C exp 




where 



C 



1 

27 



exp 



TP 3 T P 3 X(3 

3 



G 



0( 7 3 ) 
(28) 

(29) 



Integration over x yields for the absolute value £ and the 
argument £1 of L + 



£ = 

n = 



P „ (TP" 

t P 3 x(3 + x 2 ) 
6 



+ 0(7 3 ) 
0( 7 2 ) 



(30) 
(31) 



where K u is the modified Bessel function of the second 
kind of order v. Before using this result for a derivation of 
the ionization rate in the quasistatic limit, it is instructive 



to compare it to the predictions of the 'asymptotic a 
proximation' Q and the 'tunneling 
which both are based on the SPM 



ic ap- 
proximation' [6j and the 'tunneling 1 approximation' [ll| 

It is important to stress that the term proportional to 
is usually ignored, if the SPM is used. Ignoring this 



term in (|28|) and integrating over x would yield instead 
of (3H 



£spm = 



7 



\f2jrrp 



exp 



T P 



(32) 



Although the limit 7 — > of the 'asymptotic approxima- 
tion' is not easily transparent from the equations given 
in [6j , a tedious analysis yields exactly the form given in 
Eq. (f2"2")) where f2 is specified in Eq. (f3"Tj) and £ must be 
substituted with £spm from (|3"2"|) . On the other hand, to 
obtain the limit 7 — > of the 'tunneling 1 approximation' 
one should substitute £ and f2 in Eq.([22|) with 



-CTunl 
^Tunl 



7 



TqN(p 2 



exp 



in 

3 



31 



0( 7 3 



0( 7 2 ) 



(33) 
(34) 



where p — yj\ + q 2 N . 

Clearly, the 'tunneling 1 approximation' agrees with 
the 'asymptotic approximation' only for ( k 0. It is 
also evident that for ( a 1 the damping factor in the 
exponent increases with q^, whereas for the 'asymptotic 




FIG. 1: Ratio between £spm (|32|l and C (|30|) as a function of 
rp 3 /3. This ratio indicates the range of validity of the SPM. 
Since p > 1, the SPM result shows very good agreement for 
large r and starts to fail only for small values of r. The latter 
case corresponds to a small binding energy or a high intensity 
of the field. 



approximation' it remains at about — r/3. Since ( « 1 
gives the main contribution to ionization, the 'tunneling 
1 approximation' underestimates the ionization rate for 
larger kinetic momenta as is numerically proven in (l2| . 

It is also instructive to check the range of the validity of 
the SPM used to obtain the 'asymptotic approximation'. 
The ratio between £spm and £ is given by 



SPM 



/7rexp[— rp 3 /3] 



£ 



/ Tp 3 K 1/3 (Tp 3 /'S) 



(35) 



Its numerical values for different parameters is shown in 
Fig.ffl For large values of rp 3 /3 the ratio is given asymp- 
totically as 



£ 



SPM 



1 + — 



5 3 
72rp3 



f 



5F 



24«V 



(36) 



and approaches 1 in the weak-field limit. For usual laser 
parameters SPM may give an error within a few percent 
and only in the extreme case of small binding energies 
(e. g. ionization of Rydberg states) and very strong fields 
the error is significantly larger. 



D. SFA rate in the quasistatic limit 

Having obtained an exact asymptotic expression for 
\L\ 2 it is now possible to derive an analytic form of the 
SFA ionization rate Wsfa ® m tne quasistatic limit. Be- 
sides the formulas obtained in the previous two subsec- 
tions some further asymptotically exact approximations 
are, however, required. For this purpose, defining the 
azimuthal angle cj> around axis parallel to F and using 



271 



dp = dC, 



(37) 
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we rewrite Eq.([T]) as 



i 

5 p 

W SFA = / d( In (1 + q 2 N ) 2 \L\ 2 S (38) 



where 



N>N 



l*r 



(39) 



If the SPM is used, one has to use (|32j) instead of (|30|) 
in (|43p . As a consequence, one obtains a similar result, 
but function i?(r, n) has to be substituted with 



f 4r 

-Rspm(t, «) = y -j B(/9, k) exp 



2r , 



dp . (47) 



Note, the argument of the cosine in Eqs. (j23]) is pro- 
portional to 7~ 3 and leads to fast oscillations, if N and ( 
are varied. Thus the contribution of this term to the final 
result is negligibly small and it is possible to substitute 
with 2£ 2 , 



wapl 
VV SFA 



N>N 



q N (l + q 2 N ) 2 £ 2 $ . (40) 



The next step is to substitute the summation over N by 
an integral. A standard approach consists of a transfor- 
mation of the sum into an integral over qjy. This allows 
to calculate differential rates, but due to the coupling of 
<7jv and £ in p it is impossible to obtain a simple ana- 
lytical expression without the use of an expansion (e. g. 
the small kinetic momentum one, gjv <C 1)- Instead of 
the use of a double integral with respect to qpj and £ we 
rewrite (|4U|) as a double integral with respect to p and x- 
Transforming the sum into the integral over p with 



E 

N>Nq 



dp 



pr 



7(1 -c 2 )' 



(41) 



and using 



one obtains 



(1-C 2 ) 



d XP (42) 



K 5 T 



OO OO 



d PP b J d x (i + xrc'(p,xMp,x) ■ 

1 —00 

(43) 

Substitution of (f3"0")) into l|43p yields an analytical expres- 
sion for the quasistatic limit of the SFA (denoted QSFA) 



with 



W-QSFA = 10 R(t, K) 



1 2 ( T P 3 



r Ht 2 



(44) 



dp (45) 



B(p, K ) = p 8 (^) / d X (l+X 2 ) 2 Hp,x)- (46) 



Eq. is one of the central results of the present work. 
It shows that the ionization rate calculated within SFA 
is proportional to the frequency lo in the limit uj — > 0. 
Thus the SFA rate vanishes in the static limit u> = 
for all binding energies and field strengths! Clearly, this 
prediction of SFA is unphysical implying that SFA is not 
applicable in the quasistatic limit for atomic systems. 
As is also clear from the derivation, this conclusion is 
not a consequence of the usually adopted SPM or small- 
momentum approximation, since they were not adopted. 

The present finding appears to be in conflict with the 
discussion given, e. g., in where it is explicitly stressed 
that SFA (corresponding to the first-order S-matrix the- 
ory in velocity gauge) approaches the correct tunneling 
limit for lj — > (and sufficiently weak fields). However, 
in [1| (and corresponding references therein) this conclu- 
sion is reached on the basis of a derivation valid for short- 
range potentials, while the present result is obtained for 
the long ranged Coulomb potential. For short-range po- 
tentials the integral over x m (HSJ) diverges and one has 
to consider the term proportional to j 2 in (|2"5"|) . This 
term removes the divergence and the obtained limit for 
the ionization rate is now uj independent in accordance 
with the discussion in [l[ . For long-range potentials there 
is no divergence in (|4"rj|) , and thus (|4"4"|) gives the corre- 
sponding quasistatic limit of SFA for that case. Clearly, 
the agreement of the SFA rate with the one predicted by 
tunneling theories that is obtained for short-range poten- 
tials cannot be used as a measure of the validity of the 
SFA, if long-range potentials are present. However, as is 
discussed below for the specific example of hydrogen- like 
atoms, the QSFA results may be used together with tun- 
neling theories to obtain an approximate Coulomb cor- 
rection factor for SFA. 

Since for long-range potentials the SFA rate leads to 
unphysical results in the quasistatic limit, one would ex- 
pect that there is very limited interest in its explicit 
calculation. However, as is shown below, the explicit 
calculation of Wqsfa is not only useful for obtaining a 
Coulomb correction factor, but it provides also an al- 
ternative recipe for an efficient though approximate cal- 
culation of SFA rates for atomic and molecular systems 
exposed to intense laser fields in a large range of experi- 
mentally relevant laser parameters. To demonstrate this, 
calculations for hydrogen-like atoms using Eqs. (|44p to 
(|46p are discussed in the next section. 
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FIG. 2: (Color online) The proportionality coefficient R(t) — 
Wqsfa/w (|45|) is shown for the complete range of r = 
(2Eh) 3 ^ 2 /F values for the different hydrogen-like states with 
n < 3. a) For small r (corresponds to a strong field F or 
a small binding energy E^) all coefficients decrease with de- 
creasing r. To partly compensate this effect all coefficients 
R(t) are scaled by factor y/r. b) For large r (corresponds 
to a weak field F or a large binding energy E^) all coeffi- 
cients are scaled by the factor exp[— 2r/3] (in the insert the 
coefficients for m > are also shown scaled by the factor 
T~ m exp[— 2t/3]). In this limit the coefficients R(t) tend to 
those given by (T5T) . 



III. QUASISTATIC LIMIT OF SFA FOR 
HYDROGEN-LIKE ATOMS 

A. Proportionality coefficient R 

In the case of bound states of hydrogen-like atoms the 
integral over x in (HU) can be analytically calculated us- 
ing the identity 



dx 



(i + x 2 y 



(2n-3)!! 
(2n-2)!!' 



(48) 



For example, the Fourier transform of the ISo state is 
given by 



16tt 



1 



K 3/2 (1 + q%f ^ 



(49) 



Hp, x) 



7^2 



2 7 7T 



K 3 p S {l+ X 2 ) 4 



(50) 



Substitution of (|50|) into (|46[) and integration over % 
yields Bis = 1. The functions B{p) for all hydrogenic 
states with principal quantum number n < 3 are listed in 
Appendix [A] Note, for hydrogen- like states the function 
B(p) is K-independent and, therefore, the proportionality 
coefficient R is a function of r only. The evaluation of R 
according to Eq. (|45[) can then simply be performed nu- 
merically. Since the integrand is a smooth exponentially 
decaying function [this is directly evident, if the SPM ap- 
proximation is adopted as in (|47[) ] , quadrature can easily 
and very efficiently be performed with high precision. 

The proportionality coefficients R(t) for a variety of 
states of hydrogen-like atoms are shown in Fig. [2] for the 
complete range of values of the inverse field parameter r. 
In Fig. [2]a the range r < 1 is shown. For better visibility 
the function R{r)/y/r is plotted instead of R(t). It is 
worth noticing that for very small r the R values for all 
states approach 0. This is a known failure of SFA, since 
a larger field intensity to binding energy ratio (and thus 
smaller r) should clearly result in a larger and not in a 
smaller ionization rate [l3j. 

In the r range shown in Fig. [2]a in which the R values 
follow the expected behavior (decreasing R for increasing 
r), the different states vary rather differently as a func- 
tion of t. It is clearly visible that R depends mostly on 
the quantum numbers / and m and only very weakly on 
n. A different dependence is found for large values of r 
as is discussed below. 



B. Weak-field limit 

The weak- field limit F — > corresponds to t — * oo. 
Using an asymptotic expansion for the modified Bessel 
function one finds that the integrand in (145|) is propor- 
tional to exp[— 2rp 3 /3]. For large values of r the inte- 
grand decays thus rapidly as p increases. Therefore, it is 
possible to use an expansion in terms of p at p — 1. This 
procedure yields 



R nlm {r)^C^\2r)-^e-^'\ 



(51) 



The general expression for the coefficients C^f A is quite 
complicated. A very simple result occurs, however, for 
m = where C^q FA = 2n is obtained. Note, in this case 
the coefficient is ^-independent. For n < 3 the coefficients 
are given by 

^QSFA_ 9 [2n(n-Q + 2|m|-l]!! 

Cnlm ~ 2n [2n(n-Z)-l]!! ' (5 j 

Noteworthy, the I dependence is in fact limited to the 
circularity of the hydrogenic state, since I appears only 
in the form n — I. 
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FIG. 3: Ratio VKsfa/Wqsfa for the IS state of a hydrogenlike atom as a function of the frequency u) for different field strength 
F and electron binding energies Eb = k 2 /2. 



In FigJ5]b function R is shown (after multiplying it 
with e 2r / 3 to remove the exponential dependence on r) 
as a function of 1/r. The weak-field limit corresponds 
thus to 1/t — > 0. As predicted, for m = the scaled 
function R approaches C^ FA = 2n in this case. Due to 
the r' m ' factor appearing in ([51]) the high \m\ states are 
harder to ionize in the weak-field limit. It is also apparent 
from Figj2]b that the characteristic dependence on the 
quantum numbers in the weak-field limit is reached only 
for 1/r < 0.1 to 0.2. For example, down to 1/t w 0.15 
the QSFA ionization rates of the 2So and 2Po states are 
almost identical. 

It is instructive to compare the quasistatic limit of 
SFA ionization rate in weak-field limit with the well- 
known quasistatic Popov-Peremolov-Terent'ev (PPT) 
formula @, 



Wppt = \C n i\ 2 flm 



k 2 (2t) 



2n-|m|-3/2 £ ~2t/3 



(53) 



where 



\Cm\ 2 



n(n + l)l(n-l - 1)1 



fi 



lm 



(2i + \m\)\ 
2l™l|m|!(; - \m\)\ 



The ionization rates VFqsfa and Wppx both include 
the exponential term exp[— 2r/3] and the factor (2t)~I" 1 I , 
but differ in the remaining part. Introducing the ratio 



/PPT 



W PPT _ 1 2 2n - 2 \C nl \ 2 f t 



! iii 



3 p3/2 K 6n 



WqSFA W C-CJSFA y n K 5/2 p2n 

it is possible to identify four factors that prevent an 
agreement between the QSFA and PPT predictions. One 
is due to the (unphysical) ui dependence of QSFA. Also 
the constant factors that depend on the quantum num- 
bers n, I, and m differ. For example, for fixed n and 
m = QSFA predicts the same ionization rate for states 
with different I whereas the PPT formula pre dicts an I 
dependence. Then there is a constant factor (y3/7r) that 




FIG. 4: (Color online) The ratios VKsfa/VKqsfa (black solid) and W^ a /Wqsfa (red dotted) are shown as a function of the 
frequency u> for all atomic states (n < 3) and an electron binding parameter n = 0.5 a. u. at the field strength F = 0.05 a. u. 
Since both curves are in relatively good agreement with each other, the 1st step in the derivation of the QSFA discussed in 
Sec lII Dl is valid also for finite values of u). 



is, however, very close to 1. Finally, both rates differ in 
their dependence on field strength and binding energy 
which is expressed as two factors to stress the n depen- 
dence or independence. 

Since QSFA is the exact asymptotic limit of SFA, 
Eq. (|54")) can be used to derive a Coulomb-corrected SFA 
rate, Wcsfa = Qppt Wsfa- Clearly, the factor Qppt de- 
rived here explicitly for atomic hydrogen could be applied 
also to other atomic or molecular systems by performing 
the evident modifications like the introduction of effective 
quantum numbers [l(|, such as n* , I* , etc. Although the 
range of validity of VFcsfa for w ^ is not directly evi- 
dent, in contrast to M^sfa it at least reaches the tunneling 
limit. Already in the past efforts have been made to de- 
rive Coulomb-correction factors for KFR theories, but so 
far the resulting rates did not lead to convincing results 
(see [l| and references therein). Based on some approx- 
imations, A. Becker at al. [2j have proposed a Coulomb 



correction factor, C 2 . Using this factor, very good agree- 
ment is found between experimental and theoretical SFA 
ionization yields for a large number of atoms and laser 
frequencies. The comparison is, however, mostly per- 
formed on a qualitative level, since the experiments did 
not provide absolute yields and thus the theoretical and 
experimental data were adjusted at one common point. 
In addition, SFA results for atomic hydrogen (with and 
without C 2 factor) are compared to full numerical solu- 
tions of the time-dependent Schrodinger equation in Q 
and again good agreement is found (on logarithmic scale) . 

For atomic hydrogen one has C 2 = n 6n /F 2n which cor- 
responds just to the last factor in ([5^)1 . Clearly, the C 2 - 
corrected SFA rate does not approach the tunneling limit 
for uj — > 0. However, the terms missing in C 2 yield for 
lu = 0.05 a. u. and the ground state of a hydrogen atom a 
factor 0.5 - 1.2 for F = 0.05 — 0.1 a. u. This can explain 
the reasonable agreement of the C 2 -corrected SFA results 
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FIG. 5: As Fig. [5] but for field strength F = 0.1 a.u. and the electron binding parameter k = 1 a. u. (Only ratio Wsfa/Wqsfa 
is shown.) 



with the ones of ab initio calculations reported for such 
parameters in (2j. However, the deviation increases by a 
factor 10 for the CO2 laser frequency or for larger n. It 
may be noted that although [2] contains also comparisons 
with experimental data obtained with a CO2 laser, the 
present work shows that the found agreement is due to 
the fact that the comparison is made on a relative scale, 
as mentioned before. In this case the erroneous u> depen- 
dence of SFA (clearly not corrected by the C 2 factor) is, 
for example, not visible. 



C. Range of validity of QSFA 

As follows from the derivation, IFqsfa in (HU is the 
exact asymptotic form of the SFA ionization rate Wsfa in 
the limit to — ► 0. It is of course interesting to investigate 
the validity regime of QSFA for non-zero values of u>. In 
fact, as is shown now, QSFA provides for a wide range of 



parameters a good approximation to SFA even for laser 
wavelengths of around 800 nm or less. 

Fig. [3] shows the ratio Wsfa/Wqsfa for the IS state 
of a hydrogenlike atom as a function of laser frequency 
for nine different values of the inverse field parameter 
t. The variation of r is achieved by using three dif- 
ferent values for both the binding-energy related quan- 
tity k and the field intensity F. The (reference) ioniza- 
tion rate Wsfa has been calculated numerically using the 
scheme described in Sec. lII Bl All curves approach unity 
for Cij — > indicating the correctness of the derivation of 
QSFA as well as numerical consistency. In the case of 
the smallest shown value of the inverse field parameter, 
t = 0.625, one notices that the ratio shows an oscillatory 
behaviour that is due to channel closings that are not re- 
solved in QSFA. The oscillation amplitude increases with 
u>, but the ratio remains in between about 0.75 and 1.25 
in the full frequency range. Therefore, QSFA is correct 
to within 25%. If one averages over the oscillations, one 
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FIG. 6: (Color online) The exact SFA rate Wsfa, the QSFA 
rate VKqsfa, and the corrected QSFA rate Wqspa ( a U scaled 
by lo~ x ) are shown for the IS state of a hydrogenlike atom 
with k = 2 a. u. at a) F — 0.05 a. u. and b) F = 0.1a. u. 
The frequencies ui corresponding to various values of Keldysh 
parameter 7 are indicated with dashed lines. 



finds an even much better quantitative agreement be- 
tween QSFA and SFA. In view of the fact that the SFA 
rate is known to overestimate the effect of channel clos- 
ings and that these pronounced channel closing features 
mostly disappear when averaging over realistic laser pa- 
rameters (envelope, focal volume etc.), the QSFA can be 
said to provide a very accurate approximation for the 
given parameters. Note, to = 0.1a. u. corresponds to a 
laser wave length of about 450 nm and thus the shown 
frequency range covers a large range of experimentally 
relevant lasers. 

Increasing r by decreasing F (but keeping k fixed) 
leads to larger oscillation amplitudes while the oscillation 
period increases. Most importantly, the average value of 
the ratio drops with increasing w below 1. QSFA starts 
to overestimate the SFA rate. Nevertheless, the oscilla- 
tion averaged QSFA rate deviates for r = 2.5 at 800 nm 
from SFA by less than 25%. Increasing r by increas- 
ing k (for fixed F — 0.2 a. u.) decreases the oscillation 
amplitude. However, the ui averaged ratio deviates more 
from unity than for smaller binding energies. Chang- 
ing k from 0.5 a. u. (corresponding to a binding energy 
Eh w 3.5 eV) to 1.0 a. u. (13.6 eV, hydrogen atom) and 
2.0 a. u. (27.2 eV, He + ) changes the ratio at u> w 0.1a. u. 
to about 0.75 and 2.0, respectively. From the represen- 
tative examples shown in Fig. [3] one can see that these 
are general trends. A decrease of F (fixed k) leads to 
larger oscillation amplitudes and deviations of the aver- 
aged ratio from unity. This limits the applicability of 
QSFA to a smaller to range. An increase of the binding 



energy (fixed F) damps the oscillation, but increases the 
deviation from unity. Combining both results it is clear 
that QSFA works best for small binding energies and high 
field strengths and thus for small values of r. However, 
t alone is not a sufficient parameter to describe the va- 
lidity of QSFA, as can be seen from the examples shown 
for t = 2.5 and 5.0. In this example QSFA works better 
for the larger value of r that is realized by enlarging both 
E b and F. 

In Figs.0] and [5] the validity of the QSFA is investi- 
gated for different initial states of hydrogen-like atoms. 
This includes all possible states with n < 3. Fig. 2] shows 
the results for F = 0.05 a. u. and k = 0.5 a. u. These are 
the same parameters as the ones used for the IS state 
in the upper left corner of Fig. [3] The results in Fig. [5] 
were on the other hand obtained with F = 0.1 a. u. and 
k = 1.0 a. u. and correspond therefore to the ones in the 
middle of Fig. [31 Again, all ratios approach unity for 
lu — > as it should be. Comparing the results for the S 
states one notices that the oscillation amplitude increases 
with n, but the deviation of the w-averaged results is very 
similar. The same trend is visible within the P states (for 
either m = or m — 1). For a given n value the oscil- 
lations are most pronounced for I = and decrease with 
increasing I. In view of the w-averaged results the range 
of validity of the QSFA as a function of lu shows, however, 
a weaker dependence on I, but is in fact decreasing for 
increasing I. For a given n and I combination (2Po and 
2Pi, 3Po and 3Pi, or 3Do, 3Di, and 3D2) the w-averaged 
ratios indicate that the range of validity of the QSFA in- 
creases with m. One may notice the close similarity of 
the results within the series lSn, 2Pi, and 3D 2 , 2P and 
3Di, as well as 2So and 3Pi. Finally, as was the case for 
the IS state, also Figs. 0] and [5] show that a larger value 
of t decreases the oscillation amplitude and increases the 
validity regime of the QSFA. 

Fig.Hshows in addition the ratio of [see Eq. pD]) ] 

and H'qsfa- The overall good agreement with the ratio 
Wsfa/Wqsfa indicates « H^sfa- Clearly, the first 

step in deriving QSFA is well justified for finite frequen- 
cies. Especially, the highly oscillatory behavior of the 
rate due to channel closings is relatively well reproduced 
by W^p^. The main reason for the deviation between 
SFA and QSFA is thus due to step 2 of the derivation 
which smoothes out the highly oscillatory behavior of 
the rate if lu is varied. 

It is instructive to investigate the main reason for the 
failure of QSFA to reproduce SFA for large values of r. It 
turns out that for large r it is most essential to consider 
in Eq. lj2"5"l) also the terms proportional to 7 2 . This yields 
(see Appendix IB) 



TI/cor 
"QSFA 



Wq SF a (55) 

where the correction factor (for a hydrogenlike IS state) 
is given by 



C cor ( 7 ) 



exp [-f /p( 7 )] 
+ [^ 2 (7)] 6/7 ) 7/12 ' 



(56) 
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This correction significantly increases the range of va- 
lidity of the quasistatic formula. Fig. [S] shows a direct 
comparison of Wsfa and Wqg r FA (both scaled by uj^ 1 ) 
for those two case where QSFA failed most severely for 
the IS state of a hydrogenlike atom, n = 2.0 a. u. and 
F = 0.05 or 0.1a. u. Besides the oscillatory behavior of 
SFA that is also not reproduced in the corrected QSFA, 
the overall agreement is very good, even if the rate varies 
by many orders of magnitude, if u> changes from to 
0.1 a. u. The results of the corrected QSFA are not only 
interesting for improving the QSFA, but they also con- 
firm once more that the two steps made in the derivation 
of the QSFA are justified. Finally, it is worth noticing 
that the range of applicability of Eq. (j55|) is not restricted 
to a range of parameters that according to the Keldysh 
parameter 7 belongs to the quasistatic regime (7 <C 1). 
Fig. [6] shows that it works even for 7 > 1 . 



On the basis of a comparison of the QSFA result with 
the prediction of the Popov-Peremolov-Terent'ev (PPT) 
formula a Coulomb correction factor is derived. This 
factor is compared to a previously proposed one that was 
supposed to be successfully adopted in a wide range of 
calculations. It is discussed that part of this success may 
be due to the fact that the comparisons to experimental 
data was only possible on a relative scale. In this case 
a number of important terms missing in the previously 
proposed Coulomb correction factor is not visible. 

The goal of this work has been the derivation of an 
anlytical expression for the SFA in the quasistatic limit 
in the presence of long-range Coulomb interactions and a 
discussion of the resulting QSFA in comparison to SFA. 
The investigation of the validity of the SFA itself by 
comparing to the results of full solutions of the time- 
dependent Schrodinger equation is presently underway 
and will be discussed elsewhere. 



IV. CONCLUSION 

The SFA (KFR theory in velocity gauge) was studied 
both analytically and numerically in the quasistatic limit. 
The derived analytical asymptotic expression (QSFA) 
shows that in the presence of long-range Coulomb in- 
teractions and thus for ionization of neutral or positively 
charged atoms or molecules the SFA rate is proportional 
to the laser frequency in this limit. This evidently un- 
physical result indicates a break-down of the SFA. Fur- 
thermore, this result shows that in contrast to the case of 
short-range potentials the SFA rate does not converge to 
the tunneling limit for weak fields, if long-range Coulomb 
interactions are present. The analytical result is sup- 
ported by a numerical study for which an efficient scheme 
for the numerical evaluation of the SFA transition ampli- 
tude has been developed. 

Using different states of hydrogenlike atoms as an ex- 
ample, the predictions of the original SFA and the QSFA 
are compared to each other. It is found that QSFA allows 
for a rather accurate prediction of the SFA rate even for 
finite laser frequencies extending in some favorable cases 
to wavelengths of 500 nm and below. It is shown that the 
validity regime of the QSFA can even be extended using 
a correction factor that is explicitly derived for IS states. 
The large range of applicability of the QSFA is of practi- 
cal interest, since its numerical evaluation is simpler than 
the one of the original SFA rate, especially in the IR and 
far-IR frequency regime. Furthermore, it is very conve- 
nient for studies of the frequency dependence of the SFA 
rate, since the QSFA is, besides a simple proportionality 
factor, u> independent. Thus the QSFA has to be eval- 
uated for a given system and field strength only once. 
In turn, the relatively large range of laser frequencies in 
which QSFA and SFA agree demonstrates that also the 
SFA rate itself is in a wide range of laser parameters only 
proportional to to. An exception is the pronounced lu de- 
pendence due to channel closings that is not reproduced 
by QSFA. 
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APPENDIX A: FUNCTIONS B{p) FOR 
DIFFERENT STATES. 

In this Appendix function B defined by Eq. (|4"6")l is 
given explicitly for all states of hydrogen-like atoms ful- 
filling n < 3. (Note, for hydrogenlike atoms function B 
is independent of k.) 
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APPENDIX B: CORRECTION FOR LARGE r. 

In this Appendix the corrected QSFA given in Eq. (|55|) 
is derived. Introducing G = 7/3 and v = v/G, we rewrite 
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iS(u+) in ((HJ as 



iS(u+) 



/(G,x) 



where 



(Bl) 



A simple correction factor for the IS state can now 
be obtained. Indeed, for this case the main contribu- 
tion comes from /)« 1 and one can multiply VFqsfa in 



Eq. 



by 



-X+i 



f(G,x) 



3 (v + xf + l 
2i VI - G 2 v 2 



dv 



(B2) 



G cor = cx P [-(2r/3)/o( 7 )] exp[-(2r/3)/ 2 ( 7 )x 2 ] (B6) 



The real part of f(G,x) can be given using a Taylor 
expansion with respect to x as 



Re/(G, X ) = l + /o(G)+/ 2 (G)x 2 + 



(B3) 



where 



fo(G) 
hi® 



3(l + 2G ,2 )sinh~ 1 G 3Vl + G 2 



4G 3 
sinn" 1 G 
G 



AG 2 



1 



For small values of G these functions can be well approxi- 
mated by f (G) M -G 2 /10, / 2 (G) » G 2 /2 and for G < 4 
they can be fitted with good accuracy by 



/o(G) 

h{G) 



11G 2 (14 + 3G 2 ) 
~ 55(28 + 15G 2 ) + 54G 4 
5G 2 + G_ -g/s 
10 + 9G 2 



(B4) 
(B5) 



The first term in (|B6[) yields an exponential increase with 
7. The second term introduces a damping for \ > 
where d = (2r/3)/a(7). Using the approximate identity 
(valid within one percent) 



-d x 2 



(l + X 2 ) 



-dx 



71(1 

T 



e/r 



-7/12 



(B7) 



to carry out the integration over % one obtains the final 
result given in Eqs. ([55]) and ([5l 



[1] A. Becker and F. H. M. Faisal, J. Phys. B: At. Mol. Phys. 
38, Rl (2005). 

[2] A. Becker, L. Plaja, P. Moreno, M. Nurhuda, and 
F. H. M. Faisal, Phys. Rev. A 64, 023408 (2001). 

[3] T. K. Kjeldsen and L. B. Madsen, Phys. Rev. A 74, 
023407 (2006). 

[4] D. B. Milosevic, G. G. Paulus, D. Bauer, and W. Becker, 

J. Phys. B: At. Mol. Phys. 39, R203 (2006). 
[5] L. V. Keldysh, Sov. Phys. JETP 20, 1307 (1965). 
[6] H. R. Reiss, Phys. Rev. A 22, 1786 (1980). 
[7] F. H. M. Faisal, J. Phys. B: At. Mol. Phys. 6, L89 (1973). 



[8] H. R. Reiss, Prog. Quant. Electr. 16, 1 (1992). 
[9] A. M. Perelomov, V. S. Popov, and M. V. Terent'ev, 
Sov. Phys. JETP 23, 924 (1966). 
[10] M. V. Ammosov, N. B. Delone, and V. P. Krainov, 

Sov. Phys. JETP 64, 1191 (1986). 
[11] H. R. Reiss and V. P. Krainov, J. Phys. A: Math. Gen. 

36, 5575 (2003). 
[12] J. Bauer, J. Phys. A: Math. Gen. 38, 521 (2005). 
[13] J. Bauer, Phys. Rev. A 73, 023421 (2006). 



